Closing the Window on Strongly Interacting Dark Matter with IceCube 
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We use the recent results on dark matter searches of the 22-string IceCube detector to probe the 
remaining allowed window for strongly interacting dark matter in the mass range 10^ <mx <10^^ 
GeV. We calculate the expected signal in the 22-string IceCube detector from the annihilation of 
such particles captured in the Sun and compare it to the detected background. As a result, the 
remaining allowed region in the mass versus cross section parameter space is ruled out. We also 
show the expected sensitivity of the complete IceCube detector with 86 strings. 
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I. INTRODUCTION 

The search for dark matter is currently one of the 
most active fields of research in experimental astroparti- 
cle physics. Common candidates are weakly interacting 
massive particles (WIMPs), with only weak and gravi- 
tational interactions with normal matter, and which en- 
compass a variety of particle types in the mass range 
from a few tens of GeV to a few hundred TeV, where the 
upper limit is based on theoretical arguments to preserve 
unitarity pj. Among WIMP candidates are the lightest 
neutralino arising in the Minimal Supersymmetric Ex- 
tension of the Standard Model (MSSM), or the hghtest 
Kaluza-Klein mode in models of Universal Extra Dimen- 
sions (UED) 0. These candidates are thermal relics from 
the Big Bang, and since they are stable, assumed to be 
able to contribute to the dark matter content in the halos 
of galaxies. 

In a different scenario, it has been shown that super- 
massive particles can be produced in the early Universe 
and account for dark matter, independently of their in- 
teraction strength with normal matter. In order to avoid 
the mass limit imposed by unitarity constraints, one can 
model non-thermal production of super-massive parti- 
cles 0, 01 • The large mass will prevent the particle to ever 
get into thermal equilibrium with the primordial plasma. 
Super-massive dark matter candidates were coined wim- 
pzillas if assumed to interact weakly with matter and 
simpzillas [sj if the interaction is strong. The produc- 
tion mechanism discussed in [3, |3| favours particles with 
masses of the order of the inflaton mass (^10"'^^ GeV), 
but particles of masses as low as a few hundred GeV can 
also be produced, normally denoted as strongly interact- 
ing massive particles (SIMPs). Although SIMPs are in 
general taken to be in the mass range ^10** < mx ^10^ 
GeV, and simpzillas, with masses mx >10^ GeV, we will 
generically call them either simpzillas. 

Searches for simpzillas have been carried out in a va- 
riety of ways in the past. Figure [T] shows the current ex- 
clusion region in the mass versus simpzilla-nucleon cross 
section parameter space. The "Mine/Space" region (blue 



square hatched) was ruled out by many different exper- 
iments underground or space-borne 0: ^- More re- 
cently, a search for simpzillas using results from direct 
detection experiments excluded most of the remaining 
parameter space for high masses (solid yellow, labeled 
"Direct") Q- A recent analysis further constrained the 
allowed region by using the Earth heat flow to set con- 
strains on the annihilation of captured dark matter in the 
Earth's center (green striped, labelled "Earth heat") [loj . 
The blank triangular regions in Figure [T] remained un- 
probed until this work. 

Whether thermally produced or not, stable dark mat- 
ter candidates can gravitationally accumulate as a halo 
in our galaxy, becoming bound in orbits in the solar sys- 
tem. Further energy losses can occur in elastic interac- 
tions with matter in large celestial bodies, like the Sun. 
If they lose enough energy they can fall below the escape 
velocity of the object, accumulating in its center and an- 
nihilating therein |lll - [l3j |. Neutrinos of energies 0(100 
GeV) can be produced from the decays of the annihi- 
lation products. This allows for an 'indirect search' for 
dark matter with neutrino telescopes, based on a search 
for an excess neutrino flux from the direction of the Sun 
over the known atmospheric neutrino background. Pre- 
dictions of simpzilla event rates and neutrino telescopes 
sensitivity have been anticipated in [l^ . 

In this paper we explore the capabilities of IceCube 
to detect high-energy neutrinos from simpzilla annihila- 
tions in the Sun. The baseline IceCube geometry con- 
sists of 80 strings with 60 digital optical modules each, 
instrumenting 1 km'^ of ice at depths between 1450 m- 
2450 m near the geographic South Pole [lB|- Recently, 
an addition of six strings forming a denser core in the 
middle of the IceCube array has been proposed in order 
to lower the neutrino energy threshold of the detector 
to about 10 GeV, the so-called "DeepCore" array [l6|. 



^ Note that part of the exclusion region shown in ^ and repro- 
duced in [id was corrected in Q 
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FIG. 1: Excluded region at 90% C.L. in the simpzilla mass 
versus cross section parameter space. The region labeled "Di- 
rect" (solid yellow region) was excluded based on direct dark 
matter detection the "Earth heat" region(green striped) 
is excluded based on the Earth's heat flow and the blue 
square-hatched "Mine/Space" region is based on many ex- 
periments, underground and space-borne 0, @]- 



Although the IceCube geometry has been optimized to 
detect ultra-high energy (>TeV) neutrinos from poten- 
tial cosmic sources with sub-degree angular resolution, 
its angular response is still adequate at 0(100 GeV) to 
perform directional searches for dark matter. IceCube 
has recently published results on the searches for WIMPs 
arising in the MSSM p3| or the simplest UED model 
using 22 strings deployed in the 2006/2007 Antarctic sea- 
son. In this work we use this 22-string detector (hereafter 
referred to as IceCube-22) which took data during the 
2007 austral winter. We set a limit in the simpzilla mass, 
mx, versus simpzilla- nucleon cross section, Uxn, phase 
space using the fact that IceCube-22 has not detected 
any neutrino excess from the direction of the Sun. We 
also estimate the sensitivity of the full IceCube detector 
in its expected final configuration with 86 strings (here- 
after referred to as IceCube-86). 



II. NEUTRINO SPECTRUM FROM 
SIMPZILLA ANNIHILATIONS 

In order to estimate the neutrino rate from simpzilla 
annihilations in the Sun, we use the capture rate, as well 
as the full-flavour neutrino flux and energy spectrum at 
the core of the Sun as determined in [6| . We then simulate 
the neutrino propagation to the Earth, including energy 
losses and oscillation effects. 

The capture rate, Tc, depends on the mass of the sim- 
pzilla, mx, and the strength of the interaction of sim- 
pzillas with nucleons, (Txn- The capture efficiency can be 
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FIG. 2: Energy distribution of the top quarks produced in 
simpzilla annihilations at production point (dashed line) and 
before decaying (full line) . The figure shows that initial-state 
gluon radiation does not play an important role in degrading 
the energy of the tops before they decay. 



is the average nucleon mass, n^ is the number density of 
nucleons in the Sun and is the radius of the Sun. For 
large cross sections, q < 1, and simpzillas lose enough 
energy in their passage through the Sun to be efficiently 
captured. The capture rate (in s^^) is then given by 



Te = 10^^(1 + y2) 
^ (240kms-i 
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On the other hand, if q>l, capture is determined by 
the relative velocity, and only low velocity simpzillas are 
trapped. In such case, the capture rate is 
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where in the above expressions 
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and nth is the speed of the Sun in the galaxy and Vq is the 
Sun escape velocity. Note that if ctxn goes to zero, then 
the capture rate goes to zero through its dependency on 
X, as it is expected. 

The simpzilla capture rate is large enough to reach 
equilibrium with the annihilation rate, F^, in the life- 
time of the solar system, and therefore we can as- 
sume T A = 
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FIG. 3: Muon-neutrino energy distribution from top decays 
at the center of the Sun (dashed) and at the IceCube location 
(full line), after energy losses and oscillations have been taken 
into account. 
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FIG. 4: Neutrino effective area of IceCube, both the 22-string 
detector and an estimation for the planned 86-string final con- 
figuration, as a function of neutrino energy. 



The neutrino spectrum from simpzilla annihilations is 
modeled as producing a pair of quarks or gluons, which 
will fragment into high multiplicity hadronic jets. Long- 
lived light and charmed hadrons will lose energy through 
interactions in the dense solar interior before decaying, 
producing low energy neutrinos. On the other hand, due 
to its short lifetime, ~ 10~^^ s, the top quark will not 
have time to hadronize or lose energy through interac- 
tion in the medium before decaying into W6. The W 
then decays into Ivi^ with similar probability into each 
flavour I — e, /i, T. Annihilations into top quarks is there- 
fore a promising channel to produce high energetic neu- 
trinos detectable by neutrino telescopes. Bottom quarks 
are also a potential source of high energy neutrinos @. 
However in this work we only consider neutrinos pro- 
duced from the top decay chain, which makes our results 
conservative. Note that although the initially available 
energy for each top is very high (the mass mx of the sim- 
pzilla), the available phase space is shared among many 
produced particles (2.8x lO^^mx/12 tops are produced 
per annihilation [6|, where mx/12 is the mass of the sim- 
pzilla in units of 10^^ GeV), so the average energy per 
produced particle is of (TeV) . The main consequence of 
this is that initial-state gluon radiation can be neglected, 
not playing an important role in further degrading the 
energy of the tops before they decay. This is illustrated 
in Figure [U which shows the energy distribution of the 
tops at the production point (dashed line) and just be- 
fore decaying (full line), obtained with a full PYTHIA [Tl| 
simulation. 

Taking into account the above considerations, the en- 
ergy distribution of the neutrinos produced from top de- 
cays at the center of the Sun can be parametrized as 



dN -t-mw 

dEu ^ ^(E, + mO[(K + mO^ - ^\\{{E, + mw)^ - m^' 

(4) 

where mw and mt are the W boson and top quark masses 
respectively, and we have omitted a constant normaliza- 
tion factor. 

The energy dependence of the neutrino flux at the de- 
tector (Earth) will be modifled by neutrino energy losses 
in their way out of the Sun and by oscillations. In order to 
include these effects we use the publicly available Wimp- 
sim code , which we have modified so it takes the neu- 
trino spectrum from equation (jj]). We use as input the 
three fiavours with their correct relative intensity from 
the top decay chain. Figure [3] shows the injection spec- 
trum of muon neutrinos at the Sun from simpzilla anni- 
hilations according to equation (jU (dashed line), and the 
resulting spectrum at the Earth (full line). The injection 
spectrum has been truncated at 10 TeV for display pur- 
poses. As mentioned earlier, the initial energy available 
in the annihilation is shared among an enormous amount 
of annihilation products and, furthermore, the energy of 
the resulting neutrinos is further reduced by energy losses 
in the Sun. The spectrum at the detector is therefore 
not very different from those of lighter dark matter can- 
didates, like MSSM neutralinos or Kaluza-Klein modes. 
This allows us to use IceCube-22 background determina- 
tion found in both WIMP and Kaluza-Klein dark matter 
analysis [13, d]. Since IceCube is mainly sensitive to 
muon neutrinos through the detection of muon tracks, 
we only use the muon-neutrino flux at the detector for 
our calculations in what follows. 
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III. RESULTS FROM ICECUBE-22 

The number of signal events, Ng, from simpzilla anni- 
hilations in the Sun predicted in a neutrino telescope of 
effective area Ae// during an exposure time T, is given 
by 

Ns{mx, ctxn) = NfBR„-rAimx,axN)-T- J Aejf dE 

(5) 

where dN^/dE is the muon neutrino flux at the detector, 
Nt is the number of tops per annihilation, BR^ is the 
branching ratio of W decaying into neutrinos (0.326) and 
T A is the annihilation rate in the center of the Sun. The 
effective area is a measure of the efficiency of the detec- 
tor and includes the neutrino- nucleon interaction prob- 
ability, the energy loss of the produced muon from the 
interaction point to the detector and the detector trigger 
and analysis efficiency. 

In order to compare our expected signal with IceCube- 
22 observations, we use the neutrino effective area of 
IceCube-22 published in [l^, and shown in Figure SI In 
104.3 days of live time IceCube-22 has detected 13 events 
from the direction of the Sun, while the expected number 
of atmospheric neutrinos is 18.5 [l^. The expected sig- 
nal strength for the given live time is shown in Figure [5l 
We have assumed a 15% uncertainty in the expectation 
of the atmospheric neutrino flux at the relevant energies 
(Ey >50 GeV) in order to calculate the 90% CL upper 
limit on a possible signal, /x^'^, following the prescription 
in [2l|. We obtain = 3.1. A scan in the (mx, CTxn) pa- 
rameter space, selecting models that predict less events 
than /j,^'^ leads to the 90% CL exclusion region shown in 
Figure [6] (read cross-hatched area). The results disfavour 
a large region of the parameter space, and we discuss 
their significance in section fVl 



IV. SENSITIVITY OF THE COMPLETED 
ICECUBE-86 DETECTOR 

The sensitivity of the completed IceCube detector can 
be estimated from the expected number of background 
events (due to atmospheric neutrinos) and the detector 
effective area. An estimation of the expected effective 
area of IceCube-86 at trigger level has been given in [l^ . 
Previous analyses of IceCube and AMANDA show that 
the effective area is reduced between a factor of 5 to 10 
between trigger level and the final analysis level, due to 
the cuts that are applied to clean the data sets and se- 
lect high quality tracks. We have taken a conservative 
approach here and we have reduced the trigger effective 
area given in [l^ by a factor between 15 for low energies 
down to 5 at the higher energies, and used that as our 
estimate of where the effective area of IceCube-86 at a 
final analysis level might be. The curve used is shown as 
the upper curve in Figured] 

We have calculated the number of atmospheric neu- 
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FIG. 5: Expected number of events in the IceCube-22 detec- 
tor, in 104.3 days of lifetime. The rate is determined by our 
simulations and is shown for five different values of the cross 
section as labeled. 



trino events in IceCube using the measured atmospheric 
neutrino flux by AMANDA [231 . We parametrized 
the curve in Figure 10 in as dN°*"/dEdf} = 
A E-T GeV-i cm-2 s'^ sr"!, with A=0. 11407 and 
7 = 3.142. The expected number of atmospheric 
neutrino events in IceCube-86 from the direction of 
the Sun in a exposure time T is then given by 
= T • / d7V^*™/d£;dr2Ae//d£; where the 
angular integration is performed over the 0.5° that the 
Sun subtends in the sky. We have assumed no contamina- 
tion of misreconstructed atmospheric muons in the final 
sample, and therefore we use only the 3.4 atmospheric 
neutrino events from the direction of the Sun predicted 
by the integral above in 365 days live time. As in the 
previous section, we have used a 15% uncertainty on the 
overall normalization of the atmospheric neutrino flux, 
as also indicated by the AMANDA measurement, and 
calculated the sensitivity including this uncertainty. The 
90% CL exclusion region from IceCube-86 is shown in 
Figure [5] as the black line. It agrees well with the predic- 
tions made in flJl. 



V. CONCLUSIONS 

We have used recently published results from the 
22-string IceCube detector on WIMP dark matter 
searches [l3, [3 to extend the search to strongly inter- 
acting dark matter candidates in the mass range ^10* < 
mx ^10^^ GeV. The fact that the IceCube results are 
compatible with the expected atmospheric neutrino back- 
ground allows us to set restrictive limits on the mass ver- 
sus nucleon cross section parameter space. 

A note on the complementarity of direct and indirect 



5 



-20 




12 14 

log Mx (GeV) 

FIG. 6: Disfavoured region in the (mx, ctxn) parameter space 
at 90% confidence level from IceCube-22 (red cross-hatched 
area), overlayed with results from previous experiments: di- 
rect searches (solid yellow area), space-borne detectors 0] 
(blue square-hatched area) and Earth heat flow analysis [10] 
(green striped area). The black line represents the 1 year 
expected sensitivity of the completed IceCube detector. 



searches is due at this point. Generally, the total cross 
section with nucleons, CTxn, can have at least two contri- 
butions, a spin-independent (or scalar) component and 
a spin-dependent (or axial-vector) component, the rel- 
ative strength of each depending on the exact assump- 
tions made on the structure and interaction of the sim- 
pzilla. There is no detailed modelling in the literature 
of the way that simpzillas interact with matter, other 
than the interaction is assumed to be "strong", as op- 



posed to weak as in the case of WIMPs. However for 
experimental searches such distinction is important be- 
cause different experiments probe different components 
of the cross section. The results from the experiments 
shown in Figure [1] refer mainly to the spin-independent 
component of the simpzilla-nucleon cross section, due to 
the nature of the targets involved or, in the case of the 
Earth heat flow analysis, due to the fact that the capture 
of dark matter candidates in the Earth proceeds mainly 
through scattering on spinless nuclei (Fe, Si, O, Mg). On 
the other hand the IceCube results are sensitive to the 
spin-dependent cross section since the Sun is the target 
for capture, which is essentially a proton target. 

In view of the previous considerations, the current re- 
sults from IceCube provide the most restrictive limits 
so far on the spin-dependent cross section of simpzillas. 
The fact that they overlap in the (mx, (Txn) parameter 
space with previous exclusion regions from experiments 
which probe mainly the spin-independent cross section, 
practically rules out simpzillas as dark matter. It would 
require an unreasonably narrow fine tuning of the mass 
and the structure of the cross section to propose a viable 
candidate that would evade all the current limits. 

In conclusion, strongly interacting heavy relics from 
the early Universe can only account for the dark matter 
if they have masses above 10^^ GeV. As these ultra 
heavy masses are not currently favored by any model, 
the window for simpzillas as the only component of dark 
matter can be considered closed. 
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